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SUMMARY 


From measured lightning- induced voltages in aircraft elec- 
trical circuits, it has been possible to derive resistive and 
inductive transfer functions relating these voltages to the 
lightning current parameters. This derivation has been accom- 
plished with a computer program, ETCAL, which was developed for 
this purpose several years ago when this analysis was first 
attempted, under NASA Contract NAS3-12019, the results of which 
were published in NASA CR-1744 (ref. 1) in 1970. The transfer 
functions obtained by the analysis of reference 1 were observed 
to be dependent upon the lightning current waveform, indicating 
a deficiency in the analytical procedure, since the airframe is 
believed to be comprised of inactive, linear elements whose char- 
acteristics remain the same regardless of current amplitude or 
waveshape. Since lightning currents are rapidly changing and of 
short duration, the time required for diffusion of such currents 
through the aircraft skin (e.g., skin effect) is probably signi- 
ficant, although it was not factored mathematically into the 
original transfer function analysis. In the belief that skin 
effect is significant, the transfer function analysis was refined 
in a continuation of this work sponsored by NASA under Contract 
NAS3-14836. A suitable representation of the time constant of 
skin penetration was obtained from the literature and factored 
into ETCAL, with the result that transfer functions now derived 
are more consistent throughout the range of lightning stroke 
waveshapes at which comparison data derived from other measure- 
ments was available. ETCAL has been re-written and is published 
in this report. With it, measured induced voltage data can be 
utilized to determine the resistive and inductive transfer func- 
tions relating lightning currents to voltages induced in an air- 
craft electrical circuit of interest. 

At the same time, work was initiated on a completely ana- 
lytical technique to arrive at the same transfer functions. This 
is a mathematical representation of an aircraft wing and a hypo- 
thetical electrical circuit conductor inside. Some simplifying 
assumptions relating to wing geometry and lightning current flow 
are made in the first attempt. The magnetic flux linking the cir- 
cuit and flux density at the conductor location are calculated 
as fimctions of an assumed lightning current filament in the wing 
skin. The contributions from a large number of such filaments 



assvuned to comprise the wing are summed to obtain the total 
magnetic flux density and flux linking the line-integral path 
of the conductor and its airframe return. From this, the trans- 
fer inductance, M, is easily derived. The resistive transfer 
function, Rg, is calculated as a function of geometry and 
material resistivity, so it is the d-c resistance and not the 
effective resistance calculated by ETCAL. Realizing this, how- 
ever, and estimating the effective resistance from it, the 
resulting values of R and M compared well with corresponding 
values obtained via EfcAL from measured induced voltages on a 
circuit inside an F89-J aircraft wing. 

The prospects for refining WING to account for physical 
peculiarities such as resistive bonds, access doors, non-uniform 
lightning current flow and others appear promising, as well as 
the possibility of modeling other sections of airframe in this 
manner and combining them to represent a complete aircraft. The 
resulting tool would be a powerful means for aircraft designers 
to evaluate prospective electrical circuit compatibility in the 
lightning environment and eliminate potential lightning effects 
problems before they become troublesome on production aircraft. 
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INTRODUCTION 


This report describes an analysis of lightning- induced volt- 
ages in aircraft electrical circuits. The work reported herein 
builds upon former analytical work performed under NASA Contract 
NAS3-12019, "Measurements and Analysis of Lightning -Induced Volt- 
ages in Aircraft Electrical Circuits" (ref. 1). In the present 
work, progress has been made towards two broad objectives: 

1. Determination of the transfer functions relating 
measured lightning- induced voltages to lightning 
current and airframe parameters. 

2. Development of a modeling technique for calculating 
lightning- induced voltages in various aircraft elec- 
trical circuits. 

The need for a greater understanding of the transfer func- 
tions of objective (1) is apparent because these functions can 
ultimately tell us more about the various factors influencing the 
susceptibility of particular aircraft electrical circuits to 
lightning. In the previous experimental program, (ref. 1) a com- 
puter program (ETCAL) was developed which enabled the reduction 
of the resistive and inductive transfer functions, Rs and M, 
respectively, from measured oscillographic induced voltage data. 

In some ciircuits upon which measurements were made, a general 
correlation between aircraft structural characteristics and these 
transfer functions was apparent, but for other circuits, the 
meaning of them was unclear. In addition, the values of Rg and M 
were often dependent upon current waveshape - a fact which caused 
some discomfort since, in theory, transfer functions should be 
dependent upon the aircraft circuit and airframe physical char- 
acteristics only. Accordingly, further analysis work was per- 
formed uijider this contract to learn more about these factors. 

If these parameters can be understood in terms of tangible wing 
structure and wing electrical circuit characteristics, aircraft 
designers may be better able to design aircraft circuits which 
minimize lightning-induced voltages. The work accomplished toward 
this objective is presented in Section I of this report. 

All research conducted in the former program was for the pur- 
pose of measuring and understanding the possible lightning -induced 
voltages in an existing airframe. Of even greater importance is 
the problem of determining the possible levels of lightning-in- 
duced voltages in new aircraft during the design stage. If a 
reliable technique were available for this purpose, aircraft de- 
signers could evaluate the consequences of various alternate 
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structural and electrical system designs as far as lightning is 
concerned. If the results showed unacceptable induced voltage 
levels, protection could be designed into the system or a remed- 
ial design modification utilized. The second objective of this 
program, therefore, was to formulate a simple analytical model 
of a basic airframe component and evaluate its authenticity. If 
successful, such a model would show the methods needed for later 
development of a more extensive model for use in aircraft design 
evaluations . 

The main text of this report is divided into two sections. 
Section I covers the transfer function analysis and development 
of the computer program, ETCAL, to derive the effective transfer 
functions from measured induced voltage data. Section II des- 
cribes the basic theory and development of a computer program, 
WING, to calculate the transfer functions for various electrical 
circuits enclosed within an aircraft wing. An Appendix includes 
the detailed mathematical derivations of the magnetic flux for- 
mulas utilized by the computer program, WING. 
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SECTION I - TRANSFER FUNCTION ANALYSIS 


TRANSFER FUNCTIONS 


A mechanism by which lightning can affect aircraft electrical 
and avionics systems is by the generation of magnetically-induced 
and resistive voltage rises within aircraft electrical circuitry. 
These voltages may or may not be harmful to the circuits them- 
selves, as well as the avionics equipment to which these circuits 
are connected. Even if the aircraft has an electrically con- 
tinuous metallic skin, its non-cylindrical geometry will enable 
some magnetic flxix to be present within the wing and fuselage, 
even if all of the lightning current were to flow through the 
skin only. This magnetic flux will link electrical circuits with- 
in these enclosures, causing induced voltages. Similarly, the 
finite resistivity of the metallic skin will permit resistive 
voltage rises within the skin (or structure) along the path of 
lightning current flow. If an aircraft electrical circuit happens 
to employ the structure as return path, then this resistive volt- 
age enters this circuit, in series with the magnetically-induced 
voltage in the same circuit and any other (normal) steady-state 
voltages present. Capacitively coupled voltages may also be pro- 
duced in these circuits; however the essentially uniform conduct- 
ing skin of metallic aircraft keeps potential differences among 
structural elements low, thereby limiting the voltages which can 
be electrostatically coupled to interior electrical circuits. In 
practice, experimental measurements have shown magnetic and resis- 
tive components to be the most predominant. 

In the previous program (ref. 1) an expression was derived 
for the induced voltage appearing at the open circuit terminals 
of any aircraft electrical circuit, as follows: 

di, (t) 


where 


R = an effective structural resistance (ohms) 
s 

an expression for the time varying lightning 
current (amperes) 

M = an effective inductive coupling factor (henrys) . 

Equation (1) permits expression of the open circuit -induced 
voltage in terms of the lightning current and two transfer 
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functions Rg and M. Equation (1) is the result of an effort to 
reduce the complete complex wing structure and circuitry to a 
simple, two terminal Thevenin equivalent circuit consisting of 
all linear components. The Thevenin equivalent circuit param- 
eters representing each wing circuit were obtained from measure- 
ments of open circuit -induced voltages (e^^) and short circuit 
currents (isc^ terminals of the wing circuits. The rea- 

son for obtaining the Thevenin equivalent circuit is that the 
induced voltage seen by any load (appropriate avionics) connected 
to the circuit may then be calculated. 

For equation (1) to be strictly correct for all time vary- 
ing functions of iL, the transfer functions Rs and M must be 
constant for all time, or at least for the time domains charac- 
teristic of natural lightning. In other words, values of Rs and 
M derived from a set of data taken at one lightning current wave- 
shape should accurately describe the relationship between another 
lightning current waveshape and the corresponding voltage induced 
in the same circuit. In previous measurements, this was not often 
the case. For example, reference 1, tables X and XIV illustrate 
Rs and M reduced from the "fast" and "slow" lightning current 
waveshapes utilized in tests on both a position lamp and armament 
power supply circuit. (Note that R is given the subscript w to 
signify its application to a wing upon which these measurements 
were made.) These tables are reproduced here as tables I and II 
and illustrate a marked change in both R„ and M when the wave- 
shape is changed from the slow (36 x 82 jusec) waveform to the 
fast (8.2 X 14 ^sec) waveform. This data is from tests on an 
F89-J aircraft wing as described in reference 1. 

This fact has been disconcerting since it diminishes the 
validity of the R and M transfer functions in predicting the volt- 
age to be induced in a circuit by a lightning current waveshape 
other than the one being applied in a test. Yet, intuition and 
basic theory of linear systems indicates that the transfer func- 
tions should not be waveshape or time dependent. In this program, 
therefore, some effort was expended to further evaluate the ini- 
tial derivation of R and M in an effort to discover the cause of 
the variations. 

Skin Effect 


In the previous analysis of reference 1, no account had been 
taken of the "skin effect" phenomena. As an introduction to skin 
effect, consider a long cylindrical conductor as shown in figure 1, 


4 



TABLE I - CALCULATED EFFECTIVE WING RESISTANCES AND MUTUAL 
INDUCTANCES FOR VARIOUS TEST CONDITIONS 


Circuit L.050, Position Light 
Conductor 2L10E18 and Airframe 


i^ Wave Form: 

Slow Wave Form 
(36 X 82 us) 

Fast Wave Form 
(8.2 X 14 us) 

Stroke 

Location 





1 

Forward End 
of Tip Tank 

" — — ■ ' 1 
450.0 

16.3 

700.0 

4.25 

4 

Outboard 

Leading 

Edge 

57.5 

0.032 

70.0 

0.5 

5 

Trailing 
Edge of 
Aileron 

25.0 

I 

2.75 

-100.0 

1.55 

7 

Center of 

Wing 

Surface 

45.0 

1.68 

- 40.0 

0.244 

10 

Inboard 

Leading 

Edge 

40.0 

0.022 

60.0 

0.31 


(Formerly table X from ref. 1) 
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TABLE II - CALCULATED EFFECTIVE WING RESISTANCES AND MUTUAL 
INDUCTANCES FOR VARIOUS TEST CONDITIONS 

Circuit S.220, Armament Power Supply 
Conductor 2SF3886E20 and Airframe 


ij^ Wave Form: 

Slow Wave Form 
(36 X 82 lis) 

Fast Wave Form 
(8.2 X 14 us) 

Stroke 

Location 




M 

(nanohenrys) 

1 

Forward End 
of Tip Tank 

0.625 

-0.00028 

2.5 

0.0177 

4 

Outboard 

Leading 

Edge 

0.50 

-0.000224 

2.5 

0.0177 

5 

Trailing 
Edge of 
Aileron 

0.75 

-0.000336 

0.88 

0.0062 

7 

Center of 

Wing 

Surface 

0.625 

-0.00028 

1.0 

-0.00798 

10 

Inboard 

Leading 

Edge 

0.50 

-0.000224 

0.25 

0.00501 


(Formerly table XIV from ref. 1) 
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FIGURE 2 -CURRENT DENSITY IN CYLINDER WALL AT VARIOUS TIMES 
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whose length is long when compared to its diameter so that end 
effects can be neglected. If a step function of current, I, is 
injected into one end of this conductor, the total current will 
initially be confined to the outer surface. 

This condition is the result of the magnetic fields, which 
initially force all the current to flow on the exterior surface 
of the cylinder, and is the phenomenon of skin effect. As time 
passes, however, the current will diffuse inward until eventually 
a condition of equilibrium is reached at tf when the current den- 
sity throughout the cylinder wall is uniform. 

As the current diffuses into the cylinder wall, the current 
density at various times would probably appear as shown in fig- 
ure 2. The current density is initially infinite and confined 
to an infinitesimally thin shell. For a final condition, tf, the 
current density will be uniform throughout the wall thickness. 

At any time the area imder the current density curve must equal 
the total input current. The current density at the inner wall 
will increase as shown in figure 3. 



FIGURE 3 - CURRENT DENSITY AT INNER WALL 
OF CYLINDER OF FIGURE 1 

Corresponding to this time variation of the current density 
must be a time varying magnetic flixx. The magnetic flux will 
initially be external to the cylinder. As the current diffuses 
through the cylinder wall, so will a magnetic field- diffuse 
through the cylinder. 
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If a conductor, AB, is placed along the center of the cylin- 
der of figure 1, and is connected to the ungrounded end of the 
cylinder as at B, a voltage, v, will be developed between this 
conductor and the grounded end of the cylinder. This voltage is 
found by integrating the electric potential around the loop 
formed by the conductor and the cylinder wall (either the inner 
or the outer wall) . 

If the current density is assumed to be uniform over the 
length of the cylinder and one integrates aroimd the loop, ABCD, 
formed by the conductor and the inner wall, the voltage v, devel- 
oped between the conductor and the cylinder at the grounded end 
can be shown to be: 


V 


i 


J^pdje 




( 2 ) 


where 

= the current density on the inner wall at any time (t). 

p - the resistivity of the conductor wall 

I = the total length of the cylinder. 

If one evaluates the voltage developed in the loop formed 
by the conductor and the outer wall of the cylinder, the path 
encircles some magnetic flux and the equation becomes; 

■ d? (3) 

where 

J 2 = the current density at the outer wall at any time (t) . 
d0 

= the time derivative of the flux buildup in the cylin- 
der wall. 

Associated with the diffusing currents is the fact that the effec- 
tive cross-section of current path is also changing. Because the 
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same total current is passing through a varying cross-sectional 
area, the effective resistance must also be varying, due to the 
relation : 

R = (4) 
a 


where 

p = material resistivity (constant) 

I 5= length (constant) 

a = area of current flow - time varying. 

By similar reasoning, the magnetic flux appearing at the 
inside of the cylinder wall is also time variant. As a result, 
the resistive and magnetic components of any voltages induced 
between an internal conductor and the cylinder wall would be 
expected to be time variant due to skin effect in addition to 
being time variant due to the lightning current waveform and its 
differential as predicted by equation (1). Since equation (1) 
is applicable for circuits using the airframe as ground return, 
and is actually the line integral taken along the circuit con- 
ductor and the inside of the airframe as return, the coupling 
factors it sees are those apparent from the inside of the struc- 
ture. These are related in turn to the current flowing on the 
inside of the structure wall and this current is itself a time 
varying function, beginning at zero and reaching a maximum some- 
time later. 

According to Witt (ref. 2) , the time constant of the build- 
up of voltage in a solid coaxial such as the cylinder of figure 1 
will be approximately 

.2 

^0^ (seconds) (5) 


where 

6 = wall thickness in meters 
p = resistivity in ohm-meters 

jLiQ = permeability of free space (4ir x 10 henrys per meter) 

The time constant thus varies with the square of the wall 
thickness of the cylinder, inversely with the resistivity of the 
material, and linearly with its permeability. Witt*s derivation 
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was for a non-magnetic cylinder and, of course, the materials 
utilized in aircraft structures are nearly all non-magnetic. The 
derivation for an obround structure such as a wing has not been 
accomplished; however, that function is not expected to differ 
much from the cylindrical case of equation (5) . 

The point to be made by this discussion is that although the 
wing is not a cylinder it does have skin thickness and there is a 
time constant associated with the diffusion of the injected cur- 
rent throughout the wing. 

The current ix, as used in equation (1) should, therefore, be 
the current appearing at the inner skin surface. In this case, 
equation (1) relating the induced voltage to the lightning current, 
iL, should be modified to express the open circuit induced volt- 
age, eocs terms of the component of the lightning current 
appearing at the inside surface of the structure, as follows: 

d(l-e"“*^) 1, (t) 

“ dt - ' - 

where : 

e = induced voltage appearing across open circuit 
terminals 

R = the effective structural resistance 
s 

M = an effective transfer inductance between the 

lightning current and the particular electrical 
circuit 

ij^(t) = lightning current (a time-varying function) 

a = the reciprocal of the time constant of current 
penetration into the aircraft skin. 

It will be noted that this is equation (1) with the factor 

(l-e‘“‘) (7) 

added to correct the total lightning current, iL, for the portion 
of it appearing on the inside of the structure. The exponent, a, 
in the above expression is the reciprocal of the time constant of 
current penetration into the aircraft skin. 
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Since much of the resistance R 3 is caused by bonds between 
various structural elements, these resistances must be lumped 
and the behavior of lightning currents flowing through them may 
not be identical to the diffusion process described for currents 
in a solid conducting skin. Therefore, the skin effect correc- 
tion factor (eq. 7) has been omitted from the resistive term in 
equation ( 6 ) and the term Rs must therefore be defined as the 
"effective” resistance to lightning current flow, rather than 
the DC resistance. 

The factor (e utilized to account for skin effect is 
not rigorously correct, but is instead an acceptable approxima- 
tion of the skin effect diffusion process - a complex mathema- 
matical function not readily expressible in terms of the aircraft 
skin geometries of interest here. The exponential function used 
instead is more mathematically tractable and well within the 
bounds of accuracy imposed by other simplifications on these 
analytical procedures. 

At time t = 0, expression (7) goes to zero and so does 
equation ( 6 ). At later times, (7) is increasing but always less 
than 1 until t = m. Thus, the effective "amount" of i^ in equa- 
tion ( 6 ) is lessened at the earlier times and thus the factor (7) 
has the overall effect of changing the rate-of-rise of i^ in the 
expression for the magnetic component. 

The magnitude of the time constant ( 1 /ct) associated with 
current penetration of the F89-J wing would vary with location 
on the wing, and would be difficult to obtain rigorously because: 

1. The wing skin thickness tapers from the center out to 
the leading and trailing edges. 

2. The wing skin thickness tapers from the root end to 
the tip. 

3. The current injection points vary. 

4. There are discontinuities such as wheel well openings, 
antenna openings, and spacing between the wing tip and 
the wing tip fuel tank. 

For the purpose of introducing skin effect into the ETCAL analysis 
to see if more consistent values of R and M would result, an aver- 
age skin thickness was assumed for the entire wing, and calcula- 
tions were run for thicknesses, 6 , of ,080" and .125" aluminum. 
Solving equation (5) for these thicknesses resulted in penetra- 
tion time constants as follows : 
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For 


.080” 

T = 
m 


aluminum, 



6p 


( 8 ) 


47t X 10 


-7 


henrys 

meter 


. 01 meter 
cm 


.080 inches 


6(2.828 microhm - centimeters) 


.3937 inches /centimeter / (9) 


f 


T = 30 microseconds, 
m 

In a similar manner, Tm for a skin thickness of 0.125” is 
74.2 microseconds. The F89-J wing skin thicknesses vary from 
0,040” to 0,267” but the materials within which most of the 
wiring is located would appear to have 0,125” as a mean. Accord- 
ingly, the time constant, 

T = 74.2 microseconds 
m 


was used to derive a, where 


a 


1 

74.2 X 10 


( 10 ) 


for use in equation (1) . 

If it is assumed that a period of 3 time constants must 
elapse before currents are equally distributed throughout the 
skin, this will be the time required for maximum current to 
appear on the inside surface of the skin. For lightning cur- 
rents rising to crest and beginning their decay in less than 
this period the skin current penetration process will have the 
effect of slowing down the rise time of lightning current reach- 
ing the inside surface of the skin and diminishing its amplitude 
there. Thus, 0.125” skins with a time constant of 74.2 micro- 
seconds should provide greater shielding of internal electrical 
circuits against the effects of the 8 and 24 microsecond rising 
lightning currents simulated in the program, than an 0.080” skin 
whose time constant is only 30 microseconds. 

ETCAL Computer Program 

The ETCAL computer program described in reference 1 was re- 
vised to factor in skin effect, as expressed in equations (5) and 
(6). This program is presented in figure 4. A flow diagram il- 
lustrating the operation of this program is presented in fig- 
ure 5. The derivation of the ETCAL analysis is otherwise the 
same as that presented in reference 1. 
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ETCAL 


1 

g 

3 

A 

5 

6 

7 

8 

9 

10 
1 1 
12 
13 
1 A 
1 5 
16 

17 

18 

19 

20 
21 
22 
100 

107 

108 
1 lO 
1 1 5 
120 
130 
1 50 
160 
165 
170 
175 


REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 


THIS PR0GRAM CALCULATES THE RESISTIVE TR/y^SFER FUNCTI0N 
R AND INDUCTIVE TRANSFER FUNCTION M FROM INDUCED 
VOLTAGES MEASURED IN A PARTICULAR CIRCUIT IN AN AIRCRAFT 
AS A FUNCTION OF THE LIGHTNING CURRENT AMPLITUDE AND 
WAVEFORM. OUTPUTS ARE: R C0HMS)* M CHENRYS>* AND A 

TABLE 0F CALCULATED INDUCED VOLTAGES EC T) VS* TIME C T) FOR 
ALL TIME^ BASED ON R AND M. EC T> IS IN V0LTS> T IN SECONDS* 
INPUTS ARE: 

II = MAX AMPLITUDE CCREST) OF LIGHTNING CURRENT CIN AMPS) 

14 = IMAX FOR EXPONENTIAL LIGHTNING CURRENT FORMULA CIN AMPS) 
A = ALPHA EXPONENT FROM THE LIGHTNING FORMULA CIN l/SECONDS) 

B = BETA ” •• 

El = MEASURED OPEN CIRCUIT INDUCED VOLTAGE AT LIGHTNING 
CURRENT CREST TIME CIN VOLTS) 

T2 = TIME OF INDUCED VOLTAGE ZERO INTERSECT CIN SECONDS) 


••STROKE L0C:**)L 


P = DURATION OF OPEN CIRCUIT INDUCED VOLTAGE CIN SECONDS) 

D = STEPS FOR THE INDUCED VOLTAGE PLOT CIN SECONDS) 

M = DATA OSCILLOGRAM NUMBER 

L = LIGHTNING STROKE ATTACHMENT POINT CANY CONVENIENT NO.) 

C = AIRCRAFT CIRCUIT IDENTIFICATION NO . CMUST BE NUMBplRS ONLY) 
S = TIME CONSTANT OF LIGHTNING CURRENT PENETRATION INTO 
THE AIRCRAFT SKIN CIN SECONDS) 

READ I 1* I4>A,B#E1^T2*P> D*M*L>C* S 
PRINT ’’ACFT CIRCUIT NOt^’l C; "OSC NO:^’)M) 

PRINT 

LET R=E1/I 1 
PRINT 

PRINT ••R=’*5R 
LET T=T2 
GO SUB 1200 
LET M=-I*R/I3 
PRINT 

PRINT ••M='*;M 
PRINT 


180 PRINT’* T CSEC) EC T) VOLTS" 

!8! PRI NT"***** + ++*+* 

190 FOR T = 0 TO P STEP D 
200 GO SUB 1200 
210 LET E=I*R+M*I3 
220 PRINT T » E 

222 IF T>P-.5E-6 THEN 8990 
230 NEXT T 

1200 LET I=CI4/I 1 )♦! 1*CEXPC-A*T)-EXPC-B*T) ) 

1300 LET I2=C I 4/M )*I 1*C-A*EXPC-A*T)+B*EXPC -B*T) ) 

1400 LET I3=C 1/S)*I*EXPC-C 1/S)+T)+C l-EXPC-C l/S)*T) )*I2 

1600 IF T<=P THEN 1900 

1700 LET 1=0 

1800 LET 12=0 

1850 LET 13=0 

1900 RETURN 

8990 PRINT 

8991 PRINT 
9000 GO TO 100 
9100 DATA 
9999 END 


FIGURE 4 - COMPUTER PROGRAM "ETCAL" IN BASIC LANGUAGE 
FOR GE TIME SHARING COMPUTER SYSTEM 
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ETCAL 



SUBROUTINE 1200 



FIGURE 5 - FLOW DIAGRAM OF COMPUTER PROGRAM, ETCAL, TO CALCULATE R 
AND M TRANSFER FUNCTIONS FROM LIGHTNING INDUCED VOLTAGE 
MEASUREMENTS. 
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Transfer Function Analysis 


In this program, measurements were made of voltages induced 
in various circuits in a complete F89-J fighter aircraft at the 
Naval Weapons Center, China Lake, California. These measure- 
ments, most of which are reported in General Electric Report No. 
SRD-72-065, "A Test Technique for Measuring Lightning-Induced 
Voltages on Aircraft Electrical Circuits" describing the experi- 
mental phase of this program, were made with low amplitude simu- 
lated lightning currents generated by the aircraft transient 
analyzer also developed under this program. The measurements at 
China Lake included circuits in which similar m.easurements had 
been made previously under the "full scale" tests described in 
reference 1 at the GE High Voltage Laboratory. Included were the 
position lamp (circuit No. L.050) and the armament power supply 
(circuit No. S.220). The values of R^ and M derived from the 
original ETCAL analysis of reference 1 are those presented in 
tables I and II of this report. 

The revised ETCAL program of figure 4 was utilized for deri- 
vation of R^ and M from the China Lake data obtained in this pro- 
gram. At China Lake, an attempt was made to duplicate lightning 
current waveshapes utilized in the work of reference 1, however, 
the aircraft transient analyzer utilized to generate the simu- 
lated lightning current at China Lake utilized a linear resistive 
waveshaping element, whereas the full-scale lightning simulator 
used earlier utilized a non-linear resistor. While fairly simi- 
lar waveshapes can be generated, sufficient differences do exist 
to affect rate-of -rise , induced voltage levels and the resulting 
ETCAL analysis. The lightning current waveshapes applied at China 
Lake included a 24 x 72 microsecond waveform (fig. 6) and a 8 x 20 
microsecond waveform (fig. 7), These figures show the actual 
lightning current waveforms applied as well as a plot of the math- 
ematical approximation used. The mathematical expression for each 
is also given in the figures, and these expressions were utilized 
in the ETCAL analysis. 

The oscillograms of open-circuit induced voltages measured 
in the position lamp and armament power supply circuits are shown 
in figure 8. The remaining data required by the ETCAL program 
were taken from these oscillograms. Table Ilia represents the 
values of the transfer functions R^ and M calculated by the pro- 
gram from data generated by the 8 x 20 microsecond and 24 x 72 
microsecond waveforms. The values of R^ and M calculated without 
consideration of skin effect for the same circuits and somewhat 
similar test conditions from reference 1 are also presented, in 
table lllb, for comparison purposes. 
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LIGHTNING CURRENT (AMPERES) 




- ACTUAL APPLIED WAVE 

- MATHEMATICAL APPROXIMATION 

= 2.3 1, [e-00208f -e-0.073t] 

I 4 = 2.3 (I, ) 

A = 0.0208 X 10® sec"' 

B = 0.073 X 10® sec’’ 



120 140 


LIGHTNING CURRENT WAVE SHAPE 



LIGHTNING CURRENT (AMPERES) 


i 

I 



TIME -MICROSECONDS 


FIGURE 7 - 8 X 20 MICROSECOND SIMULATED LIGHTNING CURRENT WAVE SHAPE 



SLOW II waveform 

(24 X 72 /is. 666 AMPERES CREST 
SEE FIGURE 5) 


FAST j|_ WAVEFORM 
(8 X 20 /IS, 1000 AMPERES CREST 
SEE FIGURE 6) 



2V0LTS/DIV 20/iS/DIV 5 VOLTS/DIV Cqc 5 /is/DIV 

CIRCUIT L.050, POSITION LAMP, WITH STROKE TO LOCATION I AT FORWARD 
END OF TIP TANK. 



0.2 V/cm Bqc 20/6/cm 0.5 V/cm 5/ts/crr 


CIRCUIT L.050, POSITION LAMP, WITH STROKE TO LOCATION 4 AT OUTBOARD 
LEADING EDGE. 



0.2 V/cm 6 qq 20 /I sec /cm 0.5 V/cm e^Q 5/isec/cm 

CIRCUIT S220, ARMAMENT POWER SUPPLY, WITH STROKE TO LOCATION I 
AT FORWARD END OF TIP TANK. 


FIGURE 8 - OPEN CIRCUIT VOLTAGES INDUCED IN THE POSITION LAMP 
a ARMAMENT POWER SUPPLY CIRCUITS OF AN F89-J 
AIRCRAFT AT NWC CHINA LAKE. 
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TABLE III - CALCULATED EFFECTIVE WING RESISTANCES AND MUTUAL 
INDUCTANCES , WITH AND WITHOUT CONSIDERATION OF 
SKIN AFFECT 


a. Calculations by ETCAL of Figure 3 with Skin Affect 
(From Measurements at China Lake on Complete F89-J Aircraft) 


CIRCUIT 

ij^ Wave Form: 

Fast Wave Form 
(8 X 20 us) 

Slow Wave Form 
(24 X 72 „s) 

Stroke 

Location 

I^QHI 




L.050 

Position 

Lamp 

(Conductor 

2L10D18) 

1 

Forward End 
of Tip Tank 

2600 

157 

1796 

121 

4 

Outboard 

Leading 

Edge 

210 

53.6 

89.8 

17.1 

S.220 

Armament 

Power 

Supply 

(Conductor 

2SF3821G20) 

1 

Forward End 
of Tip Tank 

50.0 

5.10 

59.9 

8.98 


b. Calculations by Original ETCAL of Reference 1 Without Skin Effect 
(From Measurements at GE-HVL on Wing of an F89-J Aircraft) 


CIRCUIT 

ij^ Wave Form: 

Fast Wave Form 
(8.2 X 14 iis) 

1 

Stroke 

Location 





L.050 

Position 

Lamp 

(Conductor 

2L10D18) 

1 

Forward End 
of Tip Tank 

700.0 

4.25 

450,0 

16.3 

4 

Outboard 

Leading 

Edge 

70.0 

0.5 

57.5 

0.032 

S.220 

Armament 

; Power 

Supply 

(Conductor 

2SF3886E20) 

1 

Forward End 

2.5 

0-0177 

0.625 

-0.00028 

of Tip Tank 
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Transfer Functions 


After studying table III, it should be noted that for each 
circuit and stroke location there is much less of a percentage 
change in M resulting from the change in waveshapes than there 
was without the inclusion of skin effect in the analysis. Neither 
of the values of M calculated without skin effect is the same as 
the more consistent values calculated with skin effect. While 
the values calculated without skin effect accurately describe the 
induced voltage by equation (1) for the particular waveshape from 
which they were derived, they are inaccurate for other waveshapes. 
The degree of inaccuracy varies from case to case. If, however, 

M is reduced from data by the revised ETCAL program to include 
skin effect, it is apparent that together with it will more 

adequately describe the induced voltages generated by lightning 
current waveshapes other than the one whose data they were derived 
from, when used in equation (6). 

Consistency in and M in spite of lightning current wave- 
shape variations is expected since these are fundamentally in- 
active circuit constants based on airframe material characteris- 
tics of resistivity, permeability, and geometrical distances which 
are not, at least within the range of current amplitudes consider- 
ed here, functions of applied electrical stress. The fact that 
the program calculates a fundamentally larger value of M for the 
position lamp circuit than for the armament circuit is probably 
because the position lamp employs the airframe as a return and is 
connected via a low impedance (the bulb) to the airframe at its 
extremity, whereas the armament power supply conductor upon which 
this measurement was made is open-circuited at the pylon junction 
box. In the latter case, the magnetically induced voltage in the 
loop formed by the conductor and the airframe is probably coxintered 
somewhat by a capacitively coupled component which is effectively 
in series with the magnetic component, and of opposite polarity 
to it. This is a result frequently seen in high impedance or open- 
ended circuits (ref. 1). 

When a second skin thickness of 80 mils was used to deter- 
mine the skin effect time constant by (5) , no significant changes 
resulted in and M even though this diminishes the skin penetra- 
tion time constant, Tjji, to 30 microseconds from the 74.2 micro- 
seconds applicable to a thickness of 0.125 inch. This leads to 
the supposition that the time constant, while significant, is not 
a "sensitive" function when applied in equation (6) . In view of 
the necessity to "average" many component skin thicknesses to 
arrive at a value for derivation of Tj^^ in the first place, this 
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result is fortuitous if, in fact, it holds in most cases. Addi- 
tional experience in application of the revised ETCAL program 
will have to transpD.re for this to be confirmed, however. 

The lower resistance seen for the armament power supply cir- 
cuit (S.220) when compared to the position lamp circuit (L.050) 
although for stroke location 1 at the forward end of the tip fuel 
tank in each case, can be accounted for when one considers that 
the armament circuit does not extend the full length of the wing. 
It should be remembered that the value discussed here is the 
portion of the wing skin which serves as a return for the circuit. 

It should be noted that is not defined as the d-c resis- 
tance of the structure but is an "effective" resistance instead, 
probably arising from a very complex combination of bonding and 
material resistivities themselves. No attempt has been made, as 
yet, to derive this "effective" resistance in terms of the d-c 
equivalents of these resistivities, yet it is possible to see 
some basic relationships between values of calculated R^ and 
structural characteristics. The R^^ calculated for the position 
lamp circuit for a lightning current entering the wing tip fuel 
tank is 10 to 20 times greater than that calculated for a stroke 
entering on the wing itself. This clearly indicates that the 
joint resistance between the tank and the wing adds a component 
to the induced voltage only when the lightning current flows 
through it. The joint resistance shows up in R^ derived by ETCAL 
accordingly, only when the lightning current must flow through 
this resistance. From table Ilia, the difference between R^ de- 
rived for a stroke to the tank (2600 ^ohms) and that derived for 
a stroke to the wing tip (210 ^iohms) must be the effective resis- 
tance of the joint itself. This would say that the joint has a 
resistance, by itself, of: 

2600 - 210 2390 Mohms. 

These values, from table Ilia, were derived from measure- 
ments made at China Lake on a complete F89-J aircraft, with the 
measurements made at the circuit terminals in the cockpit. The 
earlier data from measurements made on the F89-J wing alone, from 
table Illb, result in an effective joint resistance of: 

700 - 70 - 630 ^ohms. 

Different values accrue from the corresponding slow wave data 
on the same table, but the general relationships are the same. 
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Therefore, the "joint" on the F89-J at China Lake might be sup- 
posed to have a resistance three or four times that of the same 
part of the wing assembly tested earlier at the High Voltage 
Laboratory in Pittsfield. This difference may be explainable by 
the normal variation in mechanical joint bonding resistances, 
from one assembly to another, even though of the same design. 
However, the differences also might be attributable to the added 
bonding resistance of the wing - to - fuselage joint included in the 
lightning current flow path at China Lake but obviously not in- 
cluded at the High Voltage Laboratory. This bolted joint might 
have been responsible for the additional resistance of 

apparant joint at _ apparant joint 
China Lake at HVL-Pittsf ield 

( 11 ) 

(2390 /iohms) - (630 jiohms) = 1760 jiohms. 

Note that an even larger difference between corresponding 
values of M derived from the same two data sources exists, but 
this variation can b.e largely attributed to the inclusion of skin 
effect in the analysis of China Lake data. Study of equation 
(6) shows why this is so. With the factor (1-e^'^) reducing the 
amplitude of iL on the inside skin surfaces at any time t, a 
correspondingly larger value of M must appear to provide the 
same induced voltage amplitude as would be obtained from an un- 
corrected expression, as in equation (1). ETCAL calculates 
by the same method for either case so the inclusion of skin 
effect in the China Lake analysis will not affect the deriva- 
tions . 

If the wing-to-fuselage joint in the complete aircraft does 
indeed contribute to an increased R^, then it follows that a 
lightning stroke flowing from the wing tip to fuselage, or to 
the opposite wing tip, would induce higher voltages in this and 
other circuits running from the wing to the fuselage where the 
measurements were made, as compared with measurements made on 
the same circuits at the wing root only. This was, in fact, the 
case in most instances where comparisons were made between China 
Lake data and earlier measurements on the F89-J wing alone. 
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SECTION II - A MATHEMATICAL MODELING TECHNIQUE FOR CALCULATING 
PROBABLE INDUCED VOLTAGES IN AIRCRAFT ELECTRICAL 
CIRCUITS 

BASIC THEORY 


Concurrently with the analysis of experimental data just 
described, work was begun during this program on development of 
a mathematical modeling technique to enable aircraft designers 
to calculate the probable magnitudes and waveshapes of lightning- 
induced voltages in proposed airframe and electrical circuit de- 
sign and layout configurations. Such a tool would be very useful 
in enabling aircraft designers to anticipate the magnitude of the 
transfer functions associated with particular circuits and con- 
figurations associated with new aircraft designs, thereby per- 
mitting definition of the susceptibility of connected avionics 
and electric power system components to lightning- induced volt- 
ages. The need for protective measures can then be established. 

For the purpose of this objective, a wing was chosen as the 
initial airframe to be modeled since a substantial amount of ex- 
perimental data existed (ref. 1) against which to compare the 
accuracy of the model. 

From equations (1) or (5) , it is evident that if the trans- 
fer f\anctions R and M for a particular circuit within an airframe 
are known, the voltage eQ^ in the circuit can be calculated for 
any assumed lightning current waveform, ij^Ct) . 

Due to the fast rise and decay and short time duration of 
most lightning stroke currents, nearly all of this current will 
flow in the outer skins of metallic airframes rather than through 
internal spars and ribs, etc. Assuming that this is so, and also 
that the current flows uniformly in a lineal direction through a 
structure, its skin can be represented by a very large number, n, 
of infinitely small parallel current filaments. For an aircraft 
wing, this representation would be as shown in figure 9. 

Assuming that the wing circuit pictured in figure 9 is 
"grounded" to the airframe at the outboard end of the wing, the 
voltage appearing between the inboard end of the conductor and 
the inboard airframe is equal to the line integral of voltage 
induced around a plane such as ABCD in figure 9 . 
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SKIN 

y CURRENT z 



FIGURE 9 - ELECTRICAL REPRESENTATION OF WING 


If lightning current flows on the wing skin and the skin 
has a finite resistance, then a resistive voltage drop will ap- 
pear along CD and contribute to the total voltage along path 
ABCD. If iL represents the entire lightning current, presumed 
to be equally distributed among all filaments, then the resistive 
voltage drop included along path CD is given by: 


"R = if 


( 12 ) 


In (12), Rn is the resistance of one current filament, equal 
to n times the total wing skin resistance, outboard to inboard. 

The magnetically induced voltage is, by Faraday's Law, (ref. 3) 
equal to the time rate of decrease of the total magnetic flux 
linking the circuit. This is expressed as: 


dt 


where 


e = total EMF (volts) 

= total flux (webers) 
t = time ( seconds) 


(13) 
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If the voltage around the inboard circumference of the wing 
in figure 9 is assumed to be the same at all points, the voltage 
between the circuit conductor at A and the inboard circumference 
will also be the same at any point. Thus, the magnetic flux 
lines linking all planes defined by the particular wing circuit 
conductor (AB) and any longitudinal wing intersect will be equal. 
For purposes of calculation, it is convenient to choose a plane 
along the axis of symmetry (if available) . 

The portion of the magnetic flux passing through the plane 
ABCD generated by each of the skin current filaments must be 
determined and the total summarized in order to determine the 
magnetically induced voltage, e^ as defined in equation (13). 

The magnetic flux density produced at some point, p, with 
respect to a current filament is defined by the Biot-Savart Law 
(5) as: 




dl 


(14) 


where 


I = current (amperes) 

B = magnetic flux density (webers per meter'^) 
l,r = dimensions in meters ^ 

fi = permeability of the mediiam (for air = 4^ x 10 
henrys per meter) 


The total flux, <pn, passing through a given area is equal 
to the product of the area and the component of B normal to it. 
Thus , 
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( 15 ) 



where 

<p = magnetic flux (webers) contributed by each skin 
current filament 

In the wing circuit problem, the flux density, B, must be 
determined as a function of the current ij^ and applicable geom- 
etry for each current filament along the wing skin, and inte- 
grated over the surface area defined by the plane ABCD, in accor- 
dance with equation (15). 

Since the assumption has been made, for the present, that 
the total lightning current, i^, is evenly distributed among the 
n skin current filaments, the voltage e^^ can be written: 


n=n , 

j. n=n 

7 ) 

Z. dt “ 

dt 

n=l 

n=l 


(16) 


where M is expressed in henries and iL in amperes. Since e^j^ is 
related to dij^/dt by the function, 

n=n 

- y \ ^ 

n=l 


M operates as a transfer function expressing the magnetic -induced 
voltage in terms of the skin (lightning) current. 


Computer Program 

To calculate M for various wing geometries and circuit con- 
ductor positions, a computer program (WING) in BASIC language 
for use on the General Electric time-sharing computer system has 
been written. This program also calculates the resistive trans- 
fer function Rg based on skin dimensions and resistivity. 

The computer program calculates the flux density, B, per 
ampere of lightning current, and the number of flux lines (webers) 
which link the plane ABCD from each individual skin current fila- 
ment. The program performs the calculation of expression (18) for 


28 



the flux density, B, and the expression (19) for total flux lines 
as follows: 




h - i 




^ = 27 


/p2 ■ ^2 


D + L + L log 




D=D. 


D 


- D 


-‘D=D. 


(18) 


(19) 


where 


^ = magnetic flux (webers) 

D, L, i and r = dimensions in meters, as shown in figure 10 



FIGURE 10 - DEFINITION OF DIMENSIONS OF EQUATIONS (18) and (19) 

The program divides each equation by the lightning current to 
calculate the inductive transfer function, M, and the flxix den- 
sity, B, per ampere of lightning current. 
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The derivation of equations (18) and (19) from the Biot- 
Savart Law of equation (14) in terms of the definitions of fig- 
ure 10 is presented in Appendix 1 of this report. 

The computer program WING calculates the flux linking plane 
ABCD of figure 9 per ampere of skin current from each skin cur- 
rent filament by first calculating the distances of lines AB and 
CD from the filament, and identifying these as and D 2 , respec- 
tively. Equation (19) is then solved for and for D=D 2 and 

the difference taken to achieve the flux passing between and 
D 2 , which is, in fact, flux passing through the plane ABCD, as 
shown in the cross-sectional view of figure 11. 



FIGURE 11 - CROSS SECTIONAL VIEW OF WING LEADING EDGE SHOWING 
DISTANCES UTILIZED IN CALCULATION OF MAGNETIC FLUX 
PASSING THROUGH PLANE ABCD BY EQUATION (19) 


In figure 11 it is clear that the flux passing through ABCD 
generated by the skin current filament is equal to the flux, ipi, 
surrounding the filament out to a radius D^ minus the flux, ^2 5 
surrounding the filament out to D 2 . Since all of this flux is 
generated by the same fraction of lightning current, (ij^/n), it 
follows that the fluxes and 2 /J 2 can be divided by (i^/n) to 
obtain the corresponding inductive transfer function, M in hen- 
ries, as follows: 


M 

n 


^ ^ 

(Il/h) 


( 20 ) 
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In WING is designated M8 and 02 designated M7 . The 
program proceeds to calculate the flux passing through ABCD from 
each skin current filament. This contribution is given by 

«ABCD = W - M8) (21) 

for each filament. It sxims these quantities and the total is 
equal to the inductive transfer function between the wing and 
the conductor AB with airframe return. It designates this func- 
tion M, with dimensions in henrys. M is, of course, purely a 
function of geometry and, therefore, not dependent on the ampli- 
tude of lightning current. 

The flux density, B, is indeed a function of lightning cur- 
rent amplitude as well as geometry. It is calculated in accord- 
ance with equation (18) and, since it is not identical along the 
entire length, L, of the conductor AB, the program computes it 
for a point at the middle of L (half way down the wing) . A slight 
program modification would permit it to be calculated for any 
other point along AB. 

The program WING is presented in figure 12, and a flow dia- 
gram is presented in figure 13. Two subroutines are utilized in 
the program- "Subroutine 1200" calculates M and "subrouting 1000" 
calculates B . The wing geometry and other data required are 
listed at the head of the program. The geometrical quantities are 
defined by figure 14. Normally, the number of skin conductors 
required to reach a stable solution is sufficient when J is set 
equal to 1 centimeter, for a structure the size of an aircraft 
wing. More conductors may be required for accurate solutions 
near the leading and trailing edges, and in these cases values 
of J of 0.1 or 0.01 centimeters would probably provide authentic 
solutions. When the spacing between the circuit conductor AB and 
the inside of the wing becomes of the same order of magnitude as 
the spacing between skin filament conductors, the solutions cal- 
culated by the program are not considered authentic. 

COMPARISON WITH MEASURED DATA 


The program, WING, was used to calculate the transfer induc- 
tance, M, and the flux density, B, for aircraft circuit conduc- 
tors assumed to be located at various positions on the horizontal 
axis of an F89-J wing. These calculations were made in the hopes 
that some idea of the authenticity of WING could be obtained by 
comparing the transfer functions calculated by WING with those 
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W I N G 


1 RE.v] this program CALCULATES THE THAiMSFER INDUCTANCE^ M, AND 

2 REM THE TRANSFER RESISTANCE* R* RELATING INDUCED \y0LTAGE IN 

3 REM A CIRCUIT ENCLOSED BY AN AIRCRAFT WING* TD LIGHTNING 

4 REM CURRENT THROUGH THE WING. IT PRINTS D UT THE (SELECTED) 

5 REM L0CATI0N 0F THE ENCL0SED CIRCUIT C0NDUCT0R* XI* AND THE 

6 REM TRANSFER FUNCTI0NS R AND M. IT ALSU PRINTS DUT THE FLUX 

7 REM DENSITY* B* PER AMPERE 0F LIGHTNING CURRENT* AT A P0INT 

B REM 0N THE ENCL0SED CIRCUIT 0NE HALF WAY D0 WN THE WING LENGTH 

9 REM XI IS IN CENTIMETERS 

10 REM R IS IN 0HMS 

11 REM M IS IN MICR0HENRYS 

12 REM B IS IN WEBERS PER SQUARE METER PER AMPERE 0F LIGHTNING I. 

13 REM INPUTS ARE: 

14 REM SI = AVERAGE WING SKIN THICKNESS (IN CENTIMETERS) 

15 REM Rl = THE SK I N MATERIAL RESISTIVITY (IN 0HM- C ENTI M ETEHS) 

16 REM T = THE HORIZONTAL WIDTH 0F THE TRAILING EDGE (IN CM.) 

17 REM C = THE HORIZONTAL WIDTH 0F THE MAIN WING B0X (IN CM.) 

18 REM R = LEADING EDGE RADIUS (ALSO WING THICKNESS/2) (IN CM.) 

19 REM J = DISTANCE BETWEEN ASSUMED SKIN CONDUCTORS (USUALLY 

20 REM 1* .1* OR .01 CENTIMETERS) 

21 REM XI = THE LOCATION 0F THE CONDUCTOR BEING ANALYZED (IN 

22 REM CENTIMETERS FROM THE AFT END OF MAIN WING BOX) 

23 REM L = THE LENGTH OF THE WING (IN CENTIMETERS) 

100 READ S1*R1*T*C*R 

1 0 1 READ J* X 1 * L 

105 LET B=P l*R+2+C+2*SGK(Rt2+Tt2) 

106 LET N=B/J 

107 LET PI = 3. 1 41 59265 

108 LET B = 0 

109 LET B1 = 0 

1 10 LET B2 = 0 

1 1 1 LET B3 = 0 

1 12 LET M = 0 

1 1 3 LET Ml = 0 
1 1 4 LET M3 = 0 
1 1 5 LET M4=0 

1 1 6 LET M5=0 
1 1 7 LET M6 = 0 
1 18 LET M7 = 0 

119 LET M,8-0 

120 FOR 3=(J/(2*R)) TO ( ( P 1 /2) - ( J/ ( 2*R ) ) ) STEP (J/R) 

130 IF XI <C THEN 212 

140 LET D1=S0R( ( R»2)+(X1 -C) t2-2*R*(X 1-C)*C0S(0) ) 

145 LET S = 0. 5*( D1 + (X l-O+R) 

150 LET P=SOR( ( S-Dl )♦( S-(X 1-C) )♦( S-R) /S) 

1 55 LET U=P/( S- (X 1 -C) ) 

1 60 LET X=2*ATN( U) 

180 LET Z=0+X 

190 LET D2=SQR( D1 »2+(R~(Xl-C) ) »2~2*D1*(R-(X 1-C) )*COS(Z) ) 


FIGURE 12 - COMPUTER PROGRAM 'ViNG" IN BASIC LANGUAGE 
FOR GE TIME SHARING COMPUTER SYSTEM 
(Continued next page) 


32 



195 G3 SUB 1000 

196 LET B2 = B1*SINCZ+CP 1/2) > 

200 G3 SUB 1200 

2 10 NEXT 0 

211 G0 T0 220 

212 LET D2-S0RC Rt2+Rt2-2*R*R+C0SC0) ) 

214 LET X=CPl/2) -C0/2) 

216 LET D1 =S0R( CR+(C-X1 ) > t2+D2T2-2+<R+<C-XI) )+D2+C0SCX) ) 

217 G0 T0 1 9 5 

220 F0R E=CJ/2) T0 (C-CJ/2)) STEP J 
230 LET D1 = SQRCRt 2+<X 1 -E> t2) 

240 LET D2=SQR< Rt2+CR+C-E) t2) 

245 G0 SUB 1000 

246 LET B2=-B1*< CX1-E)/D1 ) 

2 50 G0 SUB 1200 

2 60 NEXT E 

270 LET W==ATNCR/T) 

280 F0R G=C J+C0SC W) )/2 T0 T- C C J*C0 SC W) ) /2 ) STEP J*C0SCW) 

290 LET Y=G*(R/T) 

300 LET Dl=S0R(Yt2+(X l+T-G) t2) 

310 LET D2=SQRC Yt2+(C+R+T-G> t2) 

315 G0 SUB 1000 

316 LET B2=-B1*< (X 1+T-G)/D1 ) 

320 G0 SUB 1200 

330 NEXT G 

340 LET A=S1*<P1*R+2*C+2*SGRC Tt2+Rt2) ) 

350 LET R2=RI*CL/A) 

360 PRINT "X 1 ="X 1 J = ’'R2 = "R2 

365 PRINT "B="B3 
370 G0 T0 1330 

1000 LET B1 =( .002/N)*C (L/2)/C D1*S0RC CL/2) t2+Dl t2) ) )*2 
1010 RETURN 

1200 LET Ml =C .002/N)*C SORC Dl 1 2 +L t 2) +L* C L0 GC SQRC D1 1 2+L t 2 ) “L ) -L0 GC D 1 ) ) ) 

1220 LET M3=C-.002/N)*Dl 

1230 LET M4=:C .002/N)*C SORC D2t2+Lt2) ) 

1240 LET M5=C .OO2/N)*C-L)*CL0GC SQRC D2t2+Lt2)+L)-L0GC D2) ) 

1250 LET M6=C - .002/N)*D2 
1260 LET M7=M1+M3 
1270 LET M8-=M4+M5+M6 
1280 LET M-M+2*CM7-M8) 

1282 LET B3=B3+2*B2 
1290 RETURN 
1 300 DATA 
1310 DATA 
1330 END 


FIGURE 12 - Continued 
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WING 



FIGURE 13 - FLOW DIAGRAM OF COMPUTER PROGRAM, WING, TO CALCULATE R AND M TRANSFER 
FUNCTIONS FOR AN AIRCRAFT WING 
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FIGURE I3-(C0NTINUED) 
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CALCULATE D, AND Do FOR FILAMENTS ON THE TRAILING EDGE 



CALCULATE CROSS-SECTIONAL 


AREA or WING SKIN 
A=S|(2¥R +2C f 2 /t2+ r2 ) 

, I 

CALCULATE 
WING RESISTANCE 

R 2 =R| (L/A) 



RGURE 13 -(CONTINUED) 
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FIGURE 14- GEOMETRICAL DATA ITEMS REQUIRED FOR WING 
COMPUTER PROGRAM 
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reduced from measured F89-J induced voltage data by ETCAL. Due 
to the basic geometrical definitions used in WING of figure 14, 
some simplification was necessary in the F89-J dimensions for 
input to the WING program. The F89-J wing is described in refer- 
ences 1 and 5. Basically, it is 680 centimeters long and tapers 
from about 46 centimeters down to about 15 centimeters thick at 
the center spar. Chord (less ailerons and flaps) distance tapers 
from 282 centimeters at the root to 193 centimeters at the tip. 
These dimensions were approximately averaged, after inspection 
of the wing, to determine the most realistic dimension inputs to 
WING. ■ 


The data inputs describing the modeled F89-J wing were: 

^1 skin thickness) = 0.318 cm 

Rh (skin material resistivity 

for aluminiam) = 2.82 x 10 ohm-cm 

T (horizontal width of trailing 

edge) = 107 cm 

C (horizontal width of main wing 

box) = 240 cm 

R (leading edge radius) = 20 cm. 


Figure 15 is a plot of the values of transfer inductance, M, 
and flux density, B, calculated by WING for circuit conductors 
located at various positions from leading to trailing edges. 

There is some question as to whether a significant portion of the 
lightning current actually flows along the ailerons and flaps as 
it does through the skin covering the main wing box and leading 
edges. In cruise flight, with ailerons and flaps in a near 0 
position, the small gaps between ailerons and flaps may be arced 
over by the lightning current, providing a continuous path. In 
many other cases, however, this may not be so. Thus, the F89-J 
wing was modeled with and without the ailerons and flaps . When 
the program was run in the latter case, the dimension, T, was 
reduced to 1 centimeter ( to avoid a discontinuity at T=0) and 
the aft spar was assumed to carry a proportionate amount of 
lightning current. Plots of M and B are presented in figure 15 
for both cases. The transfer inductance, M, is plotted against 
the left scale which is in logarithmic increments. B, the flux 
density, is plotted against the right scale which is linear about 
a zero line shown dashed across the figure. 
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FIGURE 15- CALCULATED VALUES OF TRANSFER INDUCTANCE. M. AND FLUX 
DENSITY, B, AT CONDUCTOR LOCATIONS ALONG THE AXIS OF 
SYMMETRY OF F89-J WING MODELLED BY WING COMPUTER 
PROGRAM. 
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Figure 15 shows that the value of the transfer inductance, 

M, varies from its lowest values (apparently approaching zero) 
near the leading and trailing edges, to a maximum near (but not 
at) the center of the wing box. The maximum transfer inductance 
calculated by WING is slightly more than 0.5 microhenry at loca- 
tions around 100 centimeters from the aft main spar for this 
model, but the minimum is considerably less than this in regions 
more frequently utilized for running electrical circuits, such 
as the leading edge. 

The resistance, R, is the same for all conductor locations 
since the entire airframe is involved as the return path for the 
circuit conductor, AB. This resistance only changes when the 
wing geometry changes, as when the ailerons and flaps are omitted 
from the model. The resistance values calculated were: 

R (with ailerons and flaps included) = 7.9 microhms 

R (without ailerons and flaps) - 10.35 microhms. 

It is of interest to compare the values of M of figure 15 
with those calculated from measured data by ETCAL, as presented 
on table III. 

Table III indicates values of 53.6 and 17.1 nanohenrys 
(microhenrys x 10^) for the position lamp circuit, which runs 
along the aft spar. For most of its travel it is on the inside 
(forward) of this spar, but opposite the flap it runs along the 
outside (aft) of the spar, slightly exposing it when the flap is 
down (flaps and ailerons were set about 20*^ for the China Lake 
tests during which the measured data were taken) . The location 
of the position lamp circuit forward of the spar is believed 
within a few centimeters of the spar. This circuit employs the 
airframe as return. These values of M fall within the range of 
transfer inductances calculated by WING for conductors within a 
few centimeters of the spar, when the flaps and ailerons were 
omitted from the model. This is logical since the amplitude of 
the driving voltage and currents generated by the aircraft tran- 
sient analyzer is probably insufficient to cause sparking be- 
tween flaps and ailerons, even if a full scale lightning flash 
could do so. 

In the leading edge, circuit S.220 (armament power supply) 
listed in table III may be compared with the data presented in 
figure 15, but it must be remembered that this circuit is not 
"grounded" to the airframe at the wing tip, as is the position 
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lamp circuit, and the circuit AB utilized for derivation of the 
model (eq. 13, etc). The capacitive coupling between an un- 
grounded circuit adds a component of capacitive voltage to the 
line integral around loop ABCD in figure 9, which may be of 
opposite polarity to the inductive voltage and thus appear to 
reduce the effective transfer inductance, M, calculated by ETCAL 
for table III, since ETCAL, in its present form, does not sepa- 
rate capacitive and inductive voltage components and transfer 
functions, but considers all non-resistive voltages to be the 
result of inductive coupling only. Thus, ETCAL calculates an M 
of between 5 and 9 nanohenry s, but a circuit located in about the 
same region within the leading edge of the modeled wing will have 
in the neighborhood of 50 nanohenrys, from figure 15. 

Unfortunately, no circuits existed through the center of 
the F89-J wing upon which suitable measurements could be made 
for comparison with the WING calculations of figure 15. WING 
calculates the maximum transfer inductance to be in this region, 
as shown in table III, and it would be of interest to obtain mea- 
surements on an unshielded circuit in this location utilizing 
the airframe as return, for comparison with the calculations which 
show maximum transfer inductances of around 500 nanohenrys (0.5 
microhenry). Such a value is greater than the transfer inductance 
calculated for any circuit in this or other induced voltage mea- 
surement programs, but since most electrical wiring is installed 
in the leading or trailing ^dges of wings, no measurements have 
been made in "center” circuits. 

The values of magnetic flux density, B, per ampere of light- 
ning current are also shown in figure 15. These calculations 
are made for a point half way down the wing, since flux density 
also varies along the axis of a conductor, in a wing of finite 
length. It is noteworthy that B is zero at locations where the 
transfer inductance, M, is greatest. B, of course, describes 
flux density at a point, whereas M describes total flux through- 
out an area. Thus, there is no inconsistency in this result. 
Instead, the result is telling us that for M to be greatest, the 
conductor should be located in a region where B is lowest, or 
zero. This applies for circuits utilizing the airframe as return 
only. It is expected that circuits with independent, parallel 
return paths, such as parallel pair or twisted pair circuits, 
will experience the greatest magnetic coupling and induced volt- 
ages where the flux density, B, is greatest . This is because 
only the area between conductors of the pair is now of importance 
in equation (13) . WING would show this if it were modified to 
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calculate flux between wires of a pair instead of between a single 
conductor and the airframe. 

Practically speaking, these results lead to two fundamental 
conclusions of importance to circuit and airframe designers: 

1. To minimize magnetically induced voltages in circuits 
utilizing the airframe as return, these circuits 
should be placed in the leading and trailing edges, 
as close to the wing skin as possible. 

2. To minimize magnetically- induced voltages in cir- 
cuits with independent return, or in multi-conductor 
circuits where the voltage between any of them must 
be minimized, the conductors should be placed at a 
location near the center of the wing, where B is 
minimum. 

Figure 15 shows that B does not vary significantly within 
about 40 centimeters from the well rounded leading edge. The 
program calculates a finite value of about 22 x 10~6 webers per 
square meter per ampere of lightning current in this region, 
and calculates a negative value of 5 times this level in the 
apex of the trailing edge (of ailerons and flaps). Here, the 
solution appears less stable and the validity of nximbers less than 
5 centimeters from the apex of the trailing edge is questionable. 

The reader is reminded that all of the values of B and M 
plotted on figure 15, and the succeeding discussion are for con- 
ductor positions on the horizontal axis of symmetry only. These 
values should not be construed as valid for conductors very far 
above or below this axis. The flux density would be the same for 
equal distances above or below the axis, due to symmetry, but the 
degree of this variation is as yet unknown. WING can be extended 
to calculate M and B for all locations within the wing, and this 
should be done before WING is used in a thorough design analysis. 
It is not known if a similar result pertains, for example, when 
a conductor is placed near the top or bottom of the wing, as 
happens when the conductor is moved toward the leading or trail- 
ing edges. It is probable that B does, in fact, increase above 
and below the point on the axis of symmetry where it is zero, 
but the degree of this increase is not known. 

The resistive transfer fimctions of 7.9 and 10.35 microhms 
with and without the ailerons and flaps respectively, are much 
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lower than the values presented in table III from ETCAL reduc- 
tions of induced voltages in the F89-J, by a factor of five or 
more. This may be due to one or both of the following reasons: 

1. The resistances calculated by WING are the d-c resis- 
tances, whereas those calculated by ETCAL for table III are 
effective resistances, which reflect skin effect. In the latter 
case, the effective resistance appears higher because the light- 
ning current has not thoroughly diffused itself throughout the 
wing skin. While ETCAL has been modified to account for skin 
effect in the magnetic voltage component, it does not yet do so 
for the resistive component . 


2. The resistance calculated by WING includes only the 
skin resistance itself, and not any resistances attributable to 
the tip tank joint or the wing root joint, or any other bonding 
resistances. The table III data most appropriately comparable to 
the WING calculations is for the position light circuit (L .050) 
with lightning current delivered to the outboard leading edge. 

The table Illb data, for the wing only, is probably most compar- 
able. varies from 70 microhms for the 8.2 x 14 microsecond 
waveform to 57.5 microhms for the slower 36 x 82 microsecond wave- 
form. Thus, it appears to decrease as the waveforms become long- 
er and more time is available for diffusion of lightning current 
into the skin. If the time constant of diffusion throughout a 
0.125 inch skin is in fact near 74.2 microseconds as given by 
equation (8) , then only a small percentage of the total skin 
cross-sectional area is available for conduction by the time (at 
crest = 36 ^s) that the effective is calculated by ETCAL for 
the 36 X 82 microsecond waveform, as follows : 


At t = 36 fxs. 


effective 

'^total 


(l.^-(36/74)) 

^^‘ 1 . 628 ^ 


( 22 ) 


= (1-.614) = 0.386 


SO the effective cross-sectional area is 38,6% of the total avail 
able area. Since resistance is inversely proportional to area. 
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this means that the d-c resistance calculated by WING should be 


Rdc =0.386 (Reffective> 


If Ref f active assumed to be the 57.5 microhms calculated 

by ETCAL for table Illb, then the d-c resistance should be: 

Rj = (.386) (57.5 microhms) 
dc 

= 22.2 microhms. 

This value is very close to the values of 10,35 and 7.9 
microhms calculated by WING, when the simplifications utilized 
in the WING program are recalled. 
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CONCLUDING PI SCUSSION 


Two computer programs, ETCAL and WING have been developed to 
calculate the resistance and inductive transfer functions R and M 
from measured data or a drawing-board design model of an aircraft 
wing. The model program, WING, does not as yet account for de- 
sign details such as openings, bonding or joint resistances, non- 
metallic materials or the presence of ribs and spars, etc. Yet, 
lightning current apparently flows almost entirely in the air- 
craft skins, and the F89-J wing from which comparison test data 
was measured is completely metallic (aluminum) and of relatively 
simple design. Thus, the inductive transfer functions calculated 
by WING compare favorably with those calculated from actual in- 
duced voltage measurements by ETCAL. 

Using the transfer functions, R and M, the voltage induced 
in any wing electrical circuit for which equation (6) is valid 
can be calculated for any assumed lightning current amplitude and 
waveshape. This can be done with transfer functions generated 
from previous measured data by ETCAL or from a drawing board de- 
sign as modeled by WING. A short computer program, ECAL (not to 
be confused with ETCAL) has been prepared (fig. 16) to calculate 
e(t) for all values of t for any combination of lightning current 
parameters and transfer functions, skin thickness and resistivity. 
ECAL is programmed to accept R and M values calculated directly from 
WING or ETCAL. Thus, it must (and does) factor in skin effect in 
calculating the voltage induced in a circuit, by equation (6). 

The reader is cautioned, therefore, not to regard the geometrically 
calculated M in WING as being the effective transfer inductance 
in inducing the magnetic voltage component in the circuit. It is 
only correct when utilized in equation (6) or ECAL. This is be- 
cause skin effect prevents all of the lightning current from being 
"seen" from inside of the wing, as discussed in Section I. 

The similarity of the transfer functions derived from the 
analytical wing model with those derived by ETCAL from actual 
induced voltage measurements on the F89-J points out that the 
voltages induced in this aircraft's circuits have primarily been 
due to magnetic flux existing legitimately within the aircraft 
structure as a result of its geometry and not due to leakage 
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ECAL 


1 REM THIS PR0GRAM> ECAL> CALCULATES THE V0LTAGE ECT) INDUCED BY 

2 REM LIGHTNING IN AN AIRCRAFT ELECTRICAL CIRCUIT IF THE LIGHTNING 

3 REM CURRENT AMPLITUDE AND WAVESHAPE* AND THE TRANSFER 

A REM INDUCTANCE* M* AND THE TRANSFER RESISTANCE* R* ARE KN0WN. 

5 REM THE LIGHTNING CURRENT IS EXPRESSED IN BEWLEY'S D0UBLE 

6 REM EXP0NENTIAL F0RMULA • THE PR0GRAM PRINTS 0UT A TABLE 

7 REM 0F TIMES AND C0RRESP0NDI NG INDUCED V0LTAGE 0RDINATES* 

8 REM REPRESENTING THE TIME-VARYING INDUCED V0LTAGE* E<T). 

9 REM REtUIRED INPUTS AREx 

10 REM 14 = LIGHTNING AMPLITUDE COEFFICIENT < KIL0AMPERES) 

11 REM A = FIRST LIGHTNING FORMULA EXPONENTIAL ( 1 /MI CR0 SECONDS) 

12 REM B = SEC0ND " 

13 REM C * RECIPROCAL 0F LIGHTNING CURRENT SKIN PENETRATION 

14 REM TIME CONSTANT C 1 /MI CR0 SECONDS) 

15 REM R =: TRANSFER RESISTANCE* FROM CETCAL) OR CWING) (OHMS) 

16 REM M s TRANSFER INDUCTANCE* FROM (ETCAL) OR (WING) (HENRYS) 

17 REM P s time DURATION OF DESIRED E( T) TABLE (MICROSECONDS) 

18 REM T = TIME INCREMENTS FOR DESIRED E( T) TABLE ” 

19 REM SUGGEST THAT DATA ITEMS 1 4* A* B* AND C GO IN LINE 190 AND 

20 REM ITEMS R*M*P AND T GO IN 200 FOR EASY VARIATION. 

90 READ I4*A*B*C*R*M*P* D 

100 FOR T=0 TO P STEP D 

110 LET I 1 = 1 44c ( 1E6)*((-A*EXP(-A*T))^B*EXP(-B*T) + (A+C) + EXP((-A-C)*T)) 
120 LET 12=1 44c( 1 E6)*(B+C)*EXP( ( -B-C)*T) 

1 30 LET 13=1 1-12 

140 LET E1=R*I 4*( EXP(-A*T)-EXP(-B*T) )♦( 1-EXP( (-C)*T) ) 

150 LET E2=M*I3 
160 LET E=E1-E2 
170 PRINT T* E 
180 NEXT T 
190 DATA 
200 DATA 
210 END 


FIGURE 16 - COMPUTER PROGRAM, ECAL, TO COMPUTE LIGHTNING- 

INDUCED VOLTAGE FROM LIGHTNING CURRENT PARAMETERS 

AND TRANSFER FUNCTIONS R AND M 

s 
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flux coming through openings in the structure from outside. 

This fact is significant. It implies that plugging up, alone, 
of electrical discontinuities in the airframe may not be suf- 
ficient to reduce the succeptibility of electrical circuits 
within to acceptable levels. Treatment, instead, of individual 
circuits themselves may be required; either by modifications in 
circuit geometry, return method, or conductor shielding. The 
latter options, of course, become the only ones available when 
airframe skins are non-metallic , 

ETCAL analyses were performed only on the circuits listed 
in table III. Similar analyses would probably be successful on 
most of the other circuits measured during the China Lake tests; 
however, attempts to utilize ETCAL to derive valid transfer func- 
tions from measurements on more complex circuits in this and 
other aircraft sometimes fail. ETCAL, it will be remembered 
(ref. 1), is fundamentally a method of deriving the Thevenin 
equivalent open circuit voltage source for a two- terminal pair, 
and Thevenin 's theorem states that this can be accomplished re- 
gardless of the complexity of the circuit beyond the terminals, 
so long as it consists of a combination of active sources and 
passive, linear circuit elements. In the derivation of ETCAL, 
however, discrete data values are picked from time-varying oscil- 
lograms, and the validity of the resulting R and M values calcu- 
lated depends greatly upon selection of the correct data points, 
as described in reference 1. The analysis, for example, requires 
that the point where the measured induced voltage, eoc> crosses 
the zero line (t axis) be identified. To date, this point has 
been identified on the fundamental induced voltage trace with 
normally valid results. Some induced voltages, however, begin 
with a high frequency oscillation, usually decaying very early to 
the fundamental voltage waveform. These high frequency oscilla- 
tions also cross the t axis, however, and might be considered as 
valid data points for ETCAL. Because most of these oscillatory 
events are thought to be secondary ’’ringing" in the electrical 
circuit and therefore not a part of the true induced voltage, the 
intersects from these oscillations have correctly been overlooked 
in the ETCAL analysis. It is possible, however, that some cir- 
cuits receive, as a part of the induced voltage described by 
equation (1) , a high frequency oscillatory component or some 
other component tending to obscure proper extraction of data for 
a correct ETCAL analysis. Some characteristics of circuits in 
which this seems to occur include: 
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1. circuits a large portion of which are in exposed 
areas, such as the cockpit area; 

2. circuits with parallel branches in several different 
parts of the airframe; 

3. circuits with complex filter networks included; 

4. some circuits with non-linear elements, such as 
relay cores or motor armatures, etc. 

There are other circuits which, when tested with certain 
lightning current waveforms , do not generate induced voltages 
which cross the t axis at well defined points. This may result 
when one component of the voltage is much greater than the other 
(usually the resistive component is greater than the inductive in 
these cases) or when a substantial capacitive component must be 
added in series with the resistive and magnetic components. 

In all of these cases, it is believed that proper utiliza- 
tion of ETCAL or a subsequent improvement thereof will eventually 
enable derivation of valid transfer functions for all circuits. 
Further detailed study of these peculiar cases must be under- 
taken, however, before the ETCAL-type analytical procedure will 
become clear. 

The success of this initial attempt to calculate possible 
induced voltages, via transfer functions, in wing electrical cir- 
cuits suggests that further development can provide a computer- 
ized model capable of predicting all of the induced voltage infor- 
mation required by aircraft designers for avoidance of suscepti- 
bility or development of adequate protective measures in the 
design stage. Even in its present form, WING should be useful to 
designers in comparing the possible susceptibility of various 
wiring layout options in conventional wing designs. 

Refinements which should be made in WING include handling 

of: 

1. tapered dimensions 

2. non-metallic skin sections 

3. conductors within conduits or shields 

4. separate return circuits 
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5. access doors and other openings 

6. spars and ribs 

7. more precise geometrical description of the modeled 
wing 

8. conductors above and below the axis of symmetry 

9. conductors whose length is less than that of the 
entire wing 

10. provision for unequal lightning current distributions 
and current concentration at stroke attachment points 
on the airframe. 

Some of these refinements are simple. Others will require 
more analytical work and computer time, but all will eventually 
be possible. In a similar manner, programs can be developed for 
other elements of an aircraft and added, electrically, to WING 
to enable a complete lightning-EMC analysis to be performed on 
an aircraft design. As each refinement is developed, its results 
should, where possible, be compared with results obtained from 
measured data to assure its validity. A relatively large volume 
of induced voltage data already exists for this purpose from 
recent measurements on several aircraft utilizing the aircraft 
transient analyzer. 

The increasing use of advanced, solid-state microelectronics 
in aircraft avionics, coupled with the utilization of multiplexed 
transmission techniques and central avionics integration systems 
makes adequate protection from lightning- induced voltages a must 
if mission reliability and flight safety are to be maintained. 
Future use of solid-state electric power control techniques 
(SOSTEL) may also require that greater protection against light- 
ning-induced voltages be designed into aircraft electric power 
distribution systems than is now the case. It is hoped that the 
analytical tools described herein are made use of for this pur- 
pose. The computer programs described are available from NASA 
Aerospace Safety Research and Data Institute or General Electric 
for this purpose. The programs will be kept up-to-date by the 
addition of refinements as these are developed. 
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APPENDIX I 


1. Derivation of Equations (18) and (19) of Section II. 

For the calculation of magnetic flxix density, B, at a point, 
p, at some distance r from a current carrying wire of finite 
length, L, as shown in figure A1 below: 


y 



the Biot-Savart equation (ref, 4) gives: 


_ t£L / sinQ 
- 47T ^ ^2 


(Al) 


where I = current (amperes) 

B = flux density (webers per square meter) 

L, r, ^ = dimensions (meters) 

^ = permeability (for air - 47T x 10 ^ henry s per 
meter) 

Expressing the (sinO) in terms of the geometry of figure Al, 
equation (Al) can be re-written: 
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B = ^ 

47T 




2 2 2 2 
y) + r ii-y) + r 


dy 


+ 


Ml 

4tt 




\J iy-i)^ + ((y-^)^ + r^) 


dy 


int. 


(A2) 

2^^*^ int. 


The first integral of (A2) can be re-written as equation (A3) and 
integrated by basic integral No. 173, page 71 of reference 6 as 
follows : 


1 St . . 

1 int, 


Ml 


/.2 , .,2 , „ 2 , 3/2 


dy 


.St . ^ 
1 int 


^_mI 

4it 


MI 

477 


+ y^ - 2^y + r ) 


2r(2v-2^) 


(4je^ - 4(r^i^)) \jy^ - liy + 


2 t 2 

r + L 


(A3) 


(A4) 


0 




r (y-,e)^ + 


(A5) 


where, for the basic integral No. 173 of reference 6, 

a = 1, h = -2H- and c = (r^+i^). 

The second integral of (A2) can be re-written and solved in the 
same manner as the first integral, as follows: 


«nd . . Ml 

2 “t. = ^ 


i 


(y^ - 2ly + + r^) 


072 


= Ml 
477 


Lr 







(A6) 

(A7) 


and it is seen that integral (A5) is evaluated along the y axis 
from the bottom of the wire at 0 in figure (Al) to i , and inte- 
gral (A7) is evaluated from 1 to the top of the wire at L. This 
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integration gives B as a function of position in terms of i at r 
as follows : 


BU,r) = ^ 


Lr\/e 


+ 


2 ^ 2 
+ r 




\/(L-^)^ + r^- 


(A 8 ) 


Equation (A 8 ) is thus an expression for the fltix density, B, 
at some point at a distance from the current carrying wire. This 
is equation (18) in Section II of this report. 


To find the total amount of flux, generated by this cur- 
rent, passing through an area bounded by lines parallel to the 
wire and distances and I >2 from it, it is necessary to integrate 
B(j 0 ,r) over the area defined by these lines as shown below in fig- 
ure A2 . 


y 



FIGURE A2 


For this purpose the surface integral. 



is to be solved. Substituting (A 8 ) 


dr 

in (A9) we have 


(A9) 
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dr 


int . 1 


0 - 



+ 


g-x) 


dz dr 


(AlO) 
int . 2 


Integrating first in the y direction from 0 to L, integral 
1 of (AlO) is solved by basic integral No. 128 of page 67, refer- 
ence 6 , with 

2 

x=L, a=l, c=r and n = 1 


as follows : 


int . 



(All) 


Integral 2 of (AlO) can be re-written, 


int . 



int. 2a 


(A12) 
int . 2b 


and integral 2 is now divided into integral 2a plus integral 2b, 
as shown above. Integral 2a is solved by basic integral No. 162, 
page 70 of reference 6 with 

a = 1, b = “2L, c = (L^+r^) 


as follows : 


int. 2a = 


in 

■4’^Jdi 


^ log(2X 


- 2L + 2 \ n 




.L 


2U + + r^ 


dr (A13) 


Jo 
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Evaluating integral (A13) between limits of i = 0 and L, equa- 
tion (A13) becomes : 


int. 2a = - 2L^ + + r^) dr 


(A14) 


- 1°S(-2L + 2 / 


9 9 

+ v) dr 


Simplifying, 


int. 2a = ^ - log 2 \/?^ - log(-2L + + r^) dr (A15) 


or 


int. 2a = 


Lj I 

^ 2r 


log(-2L + 2 + r^) 


dr (A16) 


or 


int. 2a = 

\ ' u 


dr 


(A17) 


i2 ^ 2 , 

L + r - J-i 


Integral 2b of (A12) is solved by basic integral No. 164, 
page 70 of reference 6 , with 

a = 1, b = -2L, c = + r^ 


int. 2b = - 2Li + 


2 

J - 2LJL + + r^ 


log 


2.e - 2L + 


iL 


(A18) 


dr 


-*o 
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Evaluating integral 2b between limits of 0 and L gives. 


Int. 2b = [^(r) + ^ log 2r - i /l 


2^2 
+ r 


L 1 

- 7 log 


-2L + 2 


+ r^l 


(A19) 


dr 


Simplifying (A19) gives. 


int . 2b = f 
417 Jd^ 


2 
1 L 


( 1 ) 


- + r^ + 7 I log: 


2r 




2 2 / J 

-2L + 2\/L + 


dr 


(A20) 


adding integral 2a and 2b by (A17) plus (A20) now gives integral 
2 of (AlO) as follows : 


int . 2 = 


Ul [% 


L , 

“ log 




2^2 ^ . 
+ r - L-* 


dr 


(A21) 


£iZ2.+ h 

477 J D ' ' r r 
iL 


log- 




t 2 a- 2 

L + r 


- 


dr 


Now, the equation for ij) can again be written from (AlO) as 

0 = int. 1 + int. 2 

(A22) 

= (All) + (A21) 

Evaluating (All) between 0 and L, combining with (A21) and sim- 
plifying gives : 
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(A23) 


or 



(A24) 


Equation (A24) must now be integrated in the r direction 
from Dt out to ^2 achieve the entire magnetic flux passing 
through the area of interest as shown in figure A2. Equation (A24) 
can be re-written as : 


_ inf 




+ L' 


dr 


.Dr 

n f ^ 

“27T 


dr 


(A25) 


The first term of this equation is integrated by basic integral 
No, 121, page 66 of reference 6 where 

2 

a = 1, and c = L , from which, 




277 



+ V 


+ L*log 



L 



(A26) 


where r is evaluated from D-|^ to D 2 . The variable D is substituted 
in (A26) for r to become equation (19) of Section II. 
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